Multi-objective resistance-capacitance optimization algorithm: An effective multi-objective algorithm for engineering design problems

Focusing on practical engineering applications, this study introduces the Multi-Objective Resistance-Capacitance Optimization Algorithm (MORCOA), a new approach for multi-objective optimization problems. MORCOA uses the transient response behaviour of resistance-capacitance circuits to navigate complex optimization landscapes and identify global optima when faced with many competing objectives. The core approach of MORCOA combines a dynamic elimination-based crowding distance mechanism with non-dominated sorting to generate an ideal and evenly distributed Pareto front. The algorithm's effectiveness is evaluated through a structured, three-phase analysis. Initially, MORCOA is applied to five benchmark problems from the ZDT test suite, with performance assessed using various metrics and compared against state-of-the-art multi-objective optimization techniques. The study then expands to include seven problems from the DTLZ benchmark collection, further validating MORCOA's effectiveness. The final phase involves applying MORCOA to six real-world constrained engineering design problems. Notably, the optimization of a honeycomb heat sink, which is crucial in thermal management systems, is a significant part of this study. This phase uses a range of performance measures to assess MORCOA's practical application and efficiency in engineering design. The results highlight MORCOA's robustness and efficiency in both real-world engineering applications and benchmark problems, demonstrating its superior capabilities compared to existing algorithms. The effective use of MORCOA in real-world engineering design problems indicates its potential as an adaptable and powerful tool for complex multi-objective optimization tasks.

Focusing on practical engineering applications, this study introduces the Multi-Objective Resistance-Capacitance Optimization Algorithm (MORCOA), a new approach for multi-objective optimization problems.MORCOA uses the transient response behaviour of resistancecapacitance circuits to navigate complex optimization landscapes and identify global optima when faced with many competing objectives.The core approach of MORCOA combines a dynamic elimination-based crowding distance mechanism with non-dominated sorting to generate an ideal and evenly distributed Pareto front.The algorithm's effectiveness is evaluated through a structured, three-phase analysis.Initially, MORCOA is applied to five benchmark problems from the ZDT test suite, with performance assessed using various metrics and compared against stateof-the-art multi-objective optimization techniques.The study then expands to include seven problems from the DTLZ benchmark collection, further validating MORCOA's effectiveness.The final phase involves applying MORCOA to six real-world constrained engineering design problems.Notably, the optimization of a honeycomb heat sink, which is crucial in thermal management systems, is a significant part of this study.This phase uses a range of performance measures

Introduction
Sophisticated computing methods called Multi-Objective Optimization (MOO) algorithms are made to address problems with numerous, generally conflicting objectives simultaneously.The MOO addresses situations where several criteria need to be taken into account and balanced, in contrast to single-objective optimization, which seeks to identify the optimum solution with respect to a single criterion.These algorithms are essential in a wide range of industries, including engineering, logistics, finance, and healthcare.These help decision-makers negotiate difficult decision-making environments and find the optimal solutions to satisfy a variety of requirements [1,2].When trade-offs between two or more competing objectives must be determined in real-world circumstances, MOO becomes necessary.For instance, it could be necessary to balance quality and cost in product design, or the objective of environmental management might be to maximize resource use while reducing environmental effects.Due to their inability to take into consideration the complex nature of the decision-making process, traditional optimization techniques that concentrate on a single target are unable to handle these issues sufficiently [3].MOO techniques are highly influential in a number of domains: (i) With the use of these algorithms, engineers may optimize designs for a variety of factors, including performance, cost, and durability.For example, maximizing a vehicle's design to maximize both safety and fuel efficiency in the automobile industry; (ii) MOO can aid in healthcare logistics by optimizing the distribution of scarce medical resources to maximize benefits across a range of patient groups; (iii) Financial portfolio optimization uses these algorithms to maximize return on investment while minimizing risk, taking into account multiple financial products with varying risk profiles and returns; (iv) They are used to balance economic benefits with environmental protection measures, such as in waste management or water resource planning, where the goal is to maximize waste recycling and minimize disposal costs while ensuring environmental sustainability; (vi) By maximizing supply chain operations across several metrics, including cost, delivery time, and carbon footprint, they help businesses adopt more effective and sustainable supply chain procedures [4,5].
The formation of optimization algorithms to tackle difficult problems has exploded in recent years.The processes seen in nature, such as biological evolution and the behaviour of different creatures, as well as physics and chemistry occurrences and ideas from human cognition, serve as some of the inspirations for these methodologies.These include the Genetic Algorithm (GA) [6], artificial bee colony [7], Particle Swarm Optimization (PSO) [8,9], ant colony optimization [10], bat algorithm [11], Differential Evolutionary (DE) algorithm [12], simulated annealing [13], harmony search [14], and cuckoo search [15], to name a few notable examples.The grey wolf optimizer [16], whale optimizer [17,18], marine predator algorithm [19], slime mould algorithm [20], equilibrium optimizer [21], dragonfly optimizer [22], firefly optimization algorithm [23], reptile search algorithm [24], African vulture optimization algorithm [25], artificial gorilla troop optimizer [26], teaching-learning-based optimization [27], moth flame optimizer [28], political optimizer [29], snake optimizer [30], spotted hyena optimizer [31], mountain gazelle optimizer [32], golden jackal optimizer [33,34], gradient-based optimizer [35], RIME algorithm [36], INFO algorithm [37], butterfly optimization algorithm [38], among other new algorithms, have also been introduced by recent developments.Improved versions of these innovative algorithms are being created to more successfully address engineering difficulties, and they have swiftly emerged as a focus point of academic research.For instance, one creative strategy blended the benefits of swarm algorithm and evolutionary optimization for engineering applications, while another method combined the mechanics of butterfly behaviour with flower pollination processes.The majority of these algorithms were first created as single-objective optimization solutions meant to tackle problems with a single objective in mind [39,40].
However, due to their complexity, real-world issues frequently require solutions that simultaneously balance a number of objectives; these are known as MOO problems (MOOPs).These issues demand specialized algorithms that can negotiate the trade-offs between conflicting objectives.Numerous engineering domains have found themselves in need of the extensive array of MOO algorithms that have been created in recent years.Depending on the viewpoints of the decision-makers concerned, each aim in these challenges usually has varying degrees of relevance [41][42][43].The way in which MOO algorithms integrate the decision-maker's choices into the optimization process can be used to categorize them: (i) Priori methods direct the search towards a certain area of the solution space by specifying preferences prior to the optimization process.(ii) Decision-makers can fine-tune their goals and trade-offs by iteratively incorporating preferences into the optimization process through interactive methods.(ii) Posteriori approaches without stating preferences beforehand, generate a set of feasible alternatives.After the optimization process, decision-makers assess these options to determine which is best [44,45].
To address MOOPs, certain strategies have adapted the MOO algorithm into single-objective algorithms through the application of the a priori technique, which involves the integration of weights predetermined by domain experts.However, this approach encounters significant challenges, primarily related to the premature assignment of weights and the consequent difficulty in generating a well-distributed Pareto Front (PF).These hurdles significantly affect the ability of algorithms to produce a PF that encompasses a comprehensive range of optimal solutions [1,46].In contrast, the a posteriori strategy offers a more robust solution by enabling the identification of a spectrum of Pareto optimal solutions within a single computational effort, obviating the need for advance directives from experts.This method is further categorized into Multi-Objective Evolutionary Algorithm (MOEA) [47], Multi-Objective S. Ravichandran et al.Cooperative Algorithm (MOCA) [48], and Multi-Objective Swarm Algorithm (MOSA) [49].A landmark development in this field was the introduction of the Vector Evaluated Genetic Algorithm (VEGA) by Schaffer in 1984, marking the inception of evolutionary algorithms' application to multi-objective contexts [50,51].Central to the operation of most MOEAs is the principle of Non-Dominated Sorting (NDS), with the Non-Dominated Sorting Genetic Algorithm (NSGA) demonstrating this approach [52].Its successor, NSGA-II, introduced enhanced NDS processes and mechanisms aimed at diversifying solutions [53].However, NSGA-II is sometimes evaluated for its tendency towards premature convergence, exposing the loss of high-quality solutions early in the process.The Strength Pareto Evolutionary Algorithm (SPEA) and its subsequent variant, SPEA-II, improved upon the truncation methodology, fitness assignment, and density evaluation strategies of its predecessor to ensure superior convergence and solution distribution [54].Similarly, the Pareto Envelope-based Selection Algorithm (PESA) and PESA-II emphasize the prioritization of non-dominated solutions, with PESA-II adopting a region-centric selection strategy [55].Despite SPEA-II's adept search capabilities, its grid structure may impair the uniformity of solution distribution.The MOEA based on Decomposition (MOEA/D) represents another essential development, segmenting MOOPs into discrete subproblems and utilizing an aggregation strategy to leverage insights from adjacent subproblems, thus optimizing the use of computational resources [56].Nonetheless, this approach might not be as effective for discontinuous PFs, presenting limitations for complex MOOPs.The MOCAs have introduced new methodologies, including parallel processing, hybrid strategies, and coevolution, to enhance optimization processes.A notable instance is the Distributed Cooperative Coevolutionary Algorithm (DCCEA), which advances solution evolution through cooperative subpopulations and an adaptive niche strategy, complemented by a solution archive, to address decision vector segmentation and rule applications [57].
Within the domain of MOSAs, which originate inspiration from natural phenomena and behaviours, Multi-Objective Particle Swarm Optimization (MOPSO) represents a significant adaptation from conventional PSO [58].MOPSO enhances PSO by integrating multi-objective functions, incorporating mutation strategies, and employing a grid mechanism to improve the distribution of the PF.However, challenges such as suboptimal convergence rates have been noted.Following MOPSO, a variety of innovative MOSAs have been introduced, each drawing on unique aspects of natural and social behaviours to address MOOPs.Among these, the Multi-Objective Grey Wolf Optimizer (MOGWO) stands out for its emulation of grey wolves' social hierarchy and hunting strategies, augmented with an external archive and strategic selection processes for effective MOOP [59,60].The Multi-Objective Ant Lion Optimizer (MOALO) [61], multi-objective exponential distribution algorithm [62], the Multi-Objective Slime Mould Algorithm (MOSMA) [63], the Multi-Objective Equilibrium Optimizer (MOEO) [64], the Multi-Objective Arithmetic Optimization Algorithm (MOAOA) [65], the Multi-Objective Moth Flame Optimizer [66][67][68], the Multi-Objective Bat Algorithm (MOBA) [69], the Multi-Objective Cuckoo Search Algorithm (MOCSA) [70], the Multi-Objective Marine Predator Algorithm (MOMPA) [71], etc. are few examples of MOSA and each offering unique strategies to enhance solution quality and PF distribution.Further advancements have seen the development of the MOSA as follows.The Multi-Objective Multi-Verse Optimizer (MOMVO) [72], the Multi-Objective Gradient-Based Optimizer (MOGBO) [73], the Multi-Objective Teaching-Learning-Based Optimizer (MOTLBO) [74,75], the Multi-Objective JAYA (MOJAYA) algorithm [76], the Multi-Objective RAO (MORAO) algorithm [77], are some other MOO algorithms.Each of these algorithms introduces novel mechanisms for navigating the complex landscape of MOPs, with a specific focus on maintaining diversity and achieving broad coverage of the solution space.
Despite these innovations, MOSAs, like their single-objective counterparts, must overcome inherent challenges such as slow convergence rates and the propensity to converge to local optima.The ultimate aim of these algorithms is to furnish decision-makers with a well-optimized PF, addressing the critical parameters of convergence, diversity, and coverage.However, the No Free Lunch (NFL) theorem posits that no single algorithm can universally excel across all MOOPs, highlighting the situational efficacy of these optimizers and underscoring the necessity for versatile, parameter-efficient MOO algorithms capable of addressing a wide array of optimization challenges [78].This highlights the continuing need for the development of adaptive, efficient MOO algorithms that require minimal parameter tuning to cater to diverse optimization scenarios.
The Resistance-Capacitance Optimization Algorithm (RCOA) draws its unique inspiration from the dynamic response of a Resistance-Capacitance (RC) circuit to abrupt changes in voltage levels [79,80].This recent mathematical-based algorithm has gathered attention for its initial application in single-objective optimization tasks, demonstrating a promising flexibility for extension into multi-objective contexts.This research paper presents a multi-objective version of RCOA, referred to as MORCOA, which has undergone evaluation through both standard benchmark tests and practical engineering MOOPs.MORCOA not only preserves the foundational characteristics of the original RCOA but also introduces a Dynamic Elimination-Based Crowding Distance (DEBCD) mechanism aimed at enhancing the diversity of solutions [81].The developments summarized in MORCOA represent significant contributions to the field of optimization, offering a robust tool for addressing complex optimization scenarios that encompass multiple, often competing, objectives.Through the integration of DEBCD along with the NDS method [53], MORCOA stands out for its ability to generate a well-distributed PF, showcasing improved performance in attempting diverse engineering challenges.The primary contributions of this research work can be expressed as follows.
• This research introduces the MORCOA, a new algorithm designed for handling MOOP.MORCOA is based on the transient response behaviour of RC circuits, enabling it to navigate complex optimization landscapes effectively and identify global optima across multiple competing objectives.• The MORCOA integrates a DEBCD mechanism with NDS to generate an optimal and well-distributed Pareto front.This approach highlights the algorithm's capability by ensuring a balanced exploration and exploitation of the search space.• The efficacy of MORCOA is rigorously evaluated through a structured three-phase analysis involving benchmark problems from both the ZDT and DTLZ test suites.This thorough evaluation process not only benchmarks MORCOA against leading MOO algorithms but also establishes its effectiveness and efficiency in addressing diverse optimization challenges.

S. Ravichandran et al.
• A significant highlight of this study is the application of MORCOA to the optimization of honeycomb heat sink designs, a critical aspect of thermal management systems.This application demonstrates MORCOA's practical utility in real-world engineering design problems, showcasing its superior performance and adaptability in optimizing the engineering design of honeycomb heat sinks, among other engineering challenges.
The structure of the paper is outlined as follows: Section 2 provides an overview of the Resistance-Capacitance Optimization Algorithm (RCOA).Section 3 details the proposed Multi-Objective Resistance-Capacitance Optimization Algorithm (MORCOA).Section 3 also introduces the basic definitions of MOO algorithms and performance metrics.Section 4 examines the performance of MORCOA across benchmark problems and its application to real-world engineering scenarios, mostly focusing on the design optimization of honeycomb heat sinks.Section 5 concludes the paper and outlines directions for future research.

Resistance-capacitance optimization algorithm
This section explains the development of the RCOA by leveraging the characteristic transient and steady-state responses of RC circuits to voltage fluctuations [80].The RCOA algorithm is ingeniously made, drawing parallels between these circuit responses and the optimization phases of exploration and exploitation.The algorithm's foundation lies in the analogy between the transient (temporary) phase and exploration phase and between the steady-state (permanent) phase and exploitation phase.The exploration phase aims to search broadly across the solution space for promising regions, whereas the exploitation phase focuses on refining these solutions to converge towards the optimum.In RCOA, the adaptation of the RC circuit's transient response facilitates exploration, while its steady-state response underpins exploitation.Achieving a balance between these phases is crucial for the algorithm's effectiveness.
The mathematical model of RCOA incorporates the transient response v o × e − t/τ , to represent exploration and the steady-state response V s × ( 1 − e − t/τ ) , for exploitation.The progression through time t, which linearly increases from 0 to 1 across iterations, is pivotal in balancing these phases.The expression for time is defined as: where it represents the current iteration, Maxit represents the maximum number of iterations, V s denotes the source voltage to the RC circuit, v o denotes the initial voltage of the capacitor, and τ denotes the time constant of the RC circuit.The analogy extends to the capacitor's charging and discharging behaviour in the RC circuit, which is modelled over a period of 5τ when the voltage source is applied or removed.Adjusting the input waveform's time or the RC circuit's time constant alters the capacitor's voltage response, thereby modulating the algorithm's convergence behaviour.The time constant τ is crucial for this modulation and is determined by: where a set to a constant value of 5 for optimal performance, and b and r 1 are random numbers within the range [0,1].This adjustment allows for the fine-tuning of the algorithm's convergence rate, directly proportional to the time constant τ.The algorithm's implementation simulates the capacitor voltage adjustment by varying the time constant or input frequency, thus establishing a dynamic relationship between voltage (v) and time (t).This principle is vital for RCOA's development, enabling the algorithm to identify and converge on optimal solutions by simulating the RC circuit's behaviour.The position update mechanism, inspired by the circuit's transient and steady-state responses, is designed to guide the algorithm towards the best solutions.The position vectors for transient-state X ts ̅→ and steady-state X ss ̅→ are updated using the equations: The final position vector is calculated as the sum of these components: This facilitates a dynamic balance between exploration (when τ < 5τ) and exploitation (when τ ≥ 5τ), optimizing the algorithm's performance across different optimization stages.The implementation procedures of RCOA are as follows.
• Define the population size N and maximum iterations Maxit.Generate initial solutions randomly within the bounds of the decision space using Eq. ( 6).
• where X denotes the current position of the population, ub denotes the upper bound, lb denotes the lower bound, and rand denotes the uniform random number in the range of 0 and 1. • Calculate the initial objective function values for the population in the multi-dimensional search space.
• Sort and archive the best objective function values to maintain a record of optimal solutions.

S. Ravichandran et al.
• Implement effective exploration and exploitation strategies to avoid local optima and improve solution quality.
• Adjust the solution positions based on the optimization process, updating the capacitor's voltage (solution values) according to the derived equations.For each member of the population i: Update the position of the current population member according to Eq. ( 5).

End For
Update the values of t and τ for the next iteration by Eq. ( 1) and Eq. ( 2), respectively.
Re-evaluate the fitness of each population member.

Update the values of Xts
̅→ and Xss ̅→ based on the new positions X(i).

End While
Return the solution corresponding to the best fitness achieved.

Proposed multi-objective resistance-capacitance optimization algorithm
This section discusses the basic concepts of the MOO algorithm, non-dominated sorting, and the crowding distance mechanisms.It also discusses the formulation procedure of the proposed MORCOA in detail.

Basic definitions
As suggested by its name, MOO involves optimizing several objectives at the same time.The amount of goals that need to be optimized is typically used to divide this domain into segments.The MOO is a type of optimization that takes into account two or three objectives.It moves into the domain of many-objective optimization, on the other hand, when four or more objectives need to be optimized.The latter is more complicated because it has more aims, which makes it outside the scope of this study but a valuable subject to explore in this realm in the future.The objective is to identify solutions that optimize a collection of objectives simultaneously.Because the objectives frequently conflict with one another, a balance or trade-off must be made between them, which leads to complexity.As opposed to single-objective optimization, which aims to find an optimal solution, this balancing seeks to locate a set of ideal solutions, each of which reflects a unique trade-off between the objectives [82,83].In mathematical terms, a MOOP can be formulated as follows: Subject to: Where x → represents the vector of decision variables, F( x → ) is the vector function of M objectives to be minimized or maximized, d denotes the problem dimensions, g i ( x → ) are inequality constraints, with m being the number of inequality constraints, h j ( x → ) are equality constraints, with p being the number of equality constraints and lb i and ub i are the lower and upper bounds for the decision variables x i .The challenge in MOO is to identify a set of solutions that are non-dominated, meaning no other solution is better for all objectives.These solutions form the PF, illustrating the trade-offs among the objectives that decision-makers must consider.This mathematical framework supports the MOO process, providing a structured approach to navigating the complex landscape of competing objectives.
In the context of MOO, the inherent challenge stems from the diversity of objective functions that represent distinct goals.These objectives, by their very nature, possess unique scales, units, and optimization directions, rendering direct comparisons between solutions within the search space impractical, if not impossible, through conventional relational operators.To address this complexity and facilitate a structured approach to MOO, the concept of Pareto dominance was introduced.This concept serves as a cornerstone for understanding and navigating the multi-dimensional trade-offs between competing objectives.Below, the study investigates defini-S.Ravichandran et al. tions of the theory of Pareto dominance and its implications for MOO.[82]-Pareto Dominance is a principle used to compare two solutions in the context of multi-objective optimization.A solution A is said to Pareto dominates another solution B if and only if A is no worse than B in all objectives and A is strictly better than B in at least one objective, assuming that all objectives are to be minimized.Mathematically, for two solutions A and B, A dominates B if, for all objectives i,

Definition 1. Pareto Dominance
, and there exists at least one objective j where ) , assuming minimization.This concept facilitates the comparison and ranking of solutions based on their performance across multiple criteria without requiring aggregation into a single objective function.Mathematically, the vector x → dominates the vector y → if and only if: The dominant relationship among y → and x → is written as follows: x → ≻ y → .
Definition 2. Pareto Optimality [82] -A solution is considered Pareto optimal if any other solution in the search space does not dominate it.Pareto optimality implies that no other solution exists that could improve some objectives without worsening at least one other objective.Therefore, a Pareto optimal solution represents a state where an objective cannot be improved without a trade-off in terms of deteriorating performance in another objective, encapsulating the essence of compromise in MOO.A solution x → is located on the Pareto optimal front if and only if: where X denotes the decision space.
Definition 3. Pareto Optimal Set [82] -The Pareto Optimal Set, also known as the Pareto Set, comprises all solutions that are Pareto optimal within the search space.This set embodies the collection of all non-dominated solutions, highlighting the diversity of potential trade-offs among the objectives that are available.Identifying the Pareto Set is a primary goal in multi-objective optimization, as it presents the decision-maker with a spectrum of viable solutions from which to choose based on their preferences and the specific context of the problem.The Pareto set is mathematically expressed as follows.
The illustration of the Pareto optimal solution is provided in Fig. 1.Fig. 1 illustrates the concept of Pareto optimal solutions in a multi-objective optimization problem, showing both the parameter space and the objective space.The parameter space represents decision variables, while the objective space plots the objective functions f 1 and f 2 .Solutions are mapped from the parameter space to the objective space, highlighting the relationship between decision variables and objectives.
The Pareto front is the boundary of non-dominated solutions, representing the best trade-offs between objectives.Dominated solutions lie within the Pareto front, indicating inferior performance, and the dominated solutions are suboptimal because better solutions exist.This visualization clarifies the trade-offs and dominance relationships in multi-objective optimization.Fig. 1 helps to understand the nature of Pareto optimality in MOO, highlighting the trade-offs between competing objectives and the concept of dominance, as well as the relationship between the parameter space and the objective space.

Performance metrics for MOO algorithms
In the context of this study, a meticulous selection of six distinct performance evaluation metrics was employed to assess the outcomes of the research thoroughly.These metrics, namely Spacing (SP), Spread (SD), Generational Distance (GD), Hypervolume (HV), Inverted Generational Distance (IGD), and Run Time (RT), each contribute unique insights into the solution's performance, collectively offering a thorough evaluation of the optimization process's efficacy [84][85][86].Fig. 2 illustrates the definitions of all metrics.Collectively, these metrics provide a comprehensive framework for evaluating the algorithm's performance, summarizing its convergence efficiency, solution diversity, and computational feasibility.
Spacing (SP): SP examines the uniformity in the spacing between adjacent solutions along the Pareto front, aiming for a balanced distribution that avoids clustering or excessive sparsity.Smaller SP values denote more uniformly spaced solutions.SP is defined as follows.
where N o denotes the number of obtained Pareto fronts and d spacing i denotes the Euclidean distances between consecutive solutions in the objective space.
Spread (SD): SD determines the solutions' distribution range on the PF, evaluating how uniform solutions populate the objective space.A diminished SD value suggests a more evenly distributed set of solutions.The SD metric is calculated as follows.
where d f denotes the distance between the first solution and the reference point in the objective space, d l denotes the distance between the last solution and the reference point in the objective space and and it denotes the Euclidean distance calculated between consecutive solutions to assess the evenness of the distribution of solutions along the Pareto front.
Generational Distance (GD): GD evaluates the average proximity of solutions generated by the algorithm to the nearest member of the true Pareto front, acting as a convergence indicator.Lower GD values denote solutions that are nearer to the optimal Pareto front.The formula for GD is provided in Eq. (15). where is the minimum sum of absolute differences across all M objective functions between solution i and the closest solution j.In other words, it is the distance between a solution in the obtained Pareto front and the closest solution in the true Pareto front.
Hypervolume (HV): HV measures the volume encompassed by the members of a Pareto front within the objective space, serving as a dual indicator of both the convergence and the diversity of solutions within the Pareto front.A superior HV value signifies a Pareto front that spans a broader area, indicative of an enhanced spread and quality of solutions.Nadir points in multi-objective optimization represent the worst objective values found among the Pareto optimal solutions.Specifically, for a Pareto front PF consisting of solutions, the nadir point s nadir is a vector where each component is the maximum value of the corresponding objective across all Pareto optimal solutions.The nadir point s nadir is defined as: where f M (s) is the value of the m-th objective function for solution s in the Pareto front PF.
where Λ denotes the Lebesgue measure (volume), PF o is the obtained Pareto front, s⪯s ʹ ⪯s nadir means that s ʹ is bounded by the solution s and the nadir point s nadir .Inverted Generational Distance (IGD): IGD, offering a converse perspective to GD, quantifies the average distance from true Pareto front points to the nearest solution identified by the algorithm.It assesses both solution convergence and diversity, with a lower IGD value indicating superior algorithm performance.The IGD computation is as follows.

Inverted Generational Distance
Here, N t represents the count of true Pareto fronts and d IGD i is the distance between a solution in the true Pareto front and the closest solution in the obtained Pareto front.
Run Time (RT): RT assesses the algorithm's computational efficiency by measuring the time required to derive solutions.In practical applications, where time and computational resources are limited, a shorter RT is indicative of a more efficient algorithm.

Non-dominated sorting
Within the framework of the MORCOA, the procedure for updating solutions adopts a nuanced approach, in contrast to its predecessor, the RCOA, where solution updates are predicated simply on fitness value improvements.Given the multi-objective optimization context, MORCOA incorporates dominance relations to refine its solution update mechanism, which is instrumental in guiding the algorithm toward the true PF.The essence of this refined approach is embedded in the application of NDS, a method that arranges solutions based on a hierarchy of dominance, where solutions less dominated are accorded higher dominance ranks [53].In MORCOA, solutions undergo an evaluation and update process articulated through NDS, which systematically assigns ranks to solutions reflecting their level of non-dominance.The following elucidates three specific scenarios under which solutions are updated or maintained within the iterative process of MORCOA based on their relative dominance rankings.
• If a candidate solution belonging to the i th population is ranked on a higher front (F) compared to the current solution's front within the same population, the current solution is supplanted by the candidate solution.For instance, if the candidate solution is positioned on front F2 and the current solution on a lower-ranked front F3, the candidate solution is adopted as the new current solution, reflecting a direct improvement in terms of dominance ranking.• In instances where both the current and candidate solutions of the i th population are positioned on the same dominance front, a random selection mechanism with a 50 % probability is employed to decide whether to update the current solution.This scenario accommodates an element of stochasticity, ensuring diversity among solutions of equivalent dominance rank.• When the dominant front of the candidate solution is ranked lower than that of the current solution for the ith population, the current solution is preserved into the subsequent generation.This decision supports the retention of superior solutions, contributing to a more focused convergence towards the PF.
The solution update strategy articulated in MORCOA, underpinned by NDS, defines a methodical approach to navigating through the complexities of multi-objective optimization.It ensures an effective balance between exploration and exploitation within the search space, thereby enhancing the algorithm's efficacy in multi-objective settings.The mathematical representation of this strategy is given in Eq. (20).where x i (t) is the position of the i th individual at iteration t, x i (t +1) is the updated position of the i th individual at iteration t + 1, v i (t +1) is the candidate position for the i th individual at iteration t + 1, p and q are indices of the fronts, with p and q being integers, F p signifies the p th front determined through non-dominated sorting, F q signifies the q th front determined through non-dominated sorting, and rand is a random number between 0 and 1.This formulation encapsulates the strategic integration of NDS within MORCOA, facilitating a structured and efficient progression towards achieving optimal multi-objective solutions.

Dynamic-elimination-based crowding distance mechanism
Crowding Distance (CD) is a pivotal metric in MOO algorithms, crucial for ensuring uniform solution distribution across the PF [87,88].CD quantifies the density of solutions surrounding a particular solution within the non-dominated set, playing an essential role in preserving solution diversity during the optimization process.To calculate CD, solutions within the non-dominated set are first organized in ascending order based on their objective function values.The CD for a specific solution, x i , is determined by the average distance to its adjacent solutions, x i+1 and x i− 1 , effectively measuring the localized solution density.Solutions positioned at the extremities of the objective function value spectrum are assigned an infinite CD to denote their critical role in defining the solution spread.A higher CD value signifies a lower density of neighbouring solutions, indicating a broader diversity among the solutions.In addition to CD, alternative strategies for managing an external archive include the ε-Pareto archiver, grid archiver, and random archiver methods [87].The grid archiver method segments the objective space into hypercubes and achieves a balanced distribution by selectively pruning solutions within these hypercubes.Conversely, the random archiver method relies on stochastic selection for solution removal from the archive.The ε-Pareto archiver method divides the objective space into ε-defined hypercubes, retaining only the solution closest to the hypercube's origin when multiple non-dominated solutions occupy the same hypercube.The CD is widely regarded as an effective, parameter-independent technique for enhancing solution diversity within MOO algorithms, facilitating the maintenance of a fixed-size external archive.This is accomplished by sequentially removing solutions with the lowest CD to manage the archive's capacity.However, the deletion of a solution necessitates the reevaluation of CDs for the remaining solutions, as the initial removal might alter the subsequent CD ranking, potentially impacting the diversity maintained within the archive.
To mitigate this issue, the implementation of a DEBCD approach is promoted [81].DEBCD allows for the targeted update of CDs for solutions directly adjacent to the removed solution, thereby minimizing the need for widespread recalculations.Specifically, upon the removal of a solution with the smallest CD, the CDs of the immediately preceding and succeeding solutions are recalculated to reflect the new configuration.This approach ensures that only the CDs of solutions in close proximity to the removed solution are adjusted, enhancing the efficiency of the archive management process.The adoption of DEBCD significantly streamlines the process of maintaining diversity within the PS when constructing an external archive.By necessitating the recalculation of CDs for a limited number of solutions directly impacted by the removal, DEBCD optimizes the computational demands of this process, which is particularly beneficial in scenarios involving the removal of multiple solutions for archive management.This methodological refinement underscores the dynamic and computational efficiency of employing DEBCD in MOO algorithms, particularly in ensuring an even and representative distribution of solutions across the PF.
Before the elimination of the solution with the lowest CD, identified as the i th solution, the CDs for the solutions immediately adjacent to it are calculated using the following formulas: For the preceding solution: For the succeeding solution: where x i denotes the i th solution in the current population, f k (x) denotes the value of the k th objective function for solution x, d k (x i− 1 ) denote the individual distance measurement based on the k th objective for the solution immediately preceding x i , d k (x i+1 ) denotes the individual distance measurement based on the k th objective for the solution immediately succeeding x i , D(x i− 1 ) denotes the aggregate crowding distance for the solution immediately preceding x i , D(x i+1 ) denotes the aggregate crowding distance for the solution immediately succeeding x i , max ( f k ) denotes the maximum value of the k th objective function in the current population, min denotes the minimum value of the k th objective function in the current population, and m denotes the total number of objective functions.Following the removal of the i th solution due to its minimal CD, the CDs for the adjacent solutions, specifically the one preceding and the one succeeding the i th solution, are recalculated and updated as follows: For recalculating the CD of the preceding solution: For recalculating the CD of the succeeding solution: The DEBCD modifies only the CDs of solutions that are directly adjacent to the removed solution, thus minimizing the need for extensive CD recalculations across the entire solution set.This method is particularly efficient as it requires recalculating the CDs for only 2n solutions when n solutions are designated for removal to maintain the external archive, significantly optimizing the computational workload.This approach, highlighted in the visualization provided in Fig. 3, showcases the practical application of DEBCD in scenarios involving two objectives, demonstrating how CDs are dynamically adjusted to maintain an evenly distributed and diverse set of solutions within the PS, thereby enhancing the efficiency of MOO algorithms.
The pseudocode outlines the procedure for managing an external archive using the DEBCD method shown in Algorithm 2, ensuring that the archive remains within a predetermined size while maintaining solution diversity through dynamic CD updates.Update the CD of the (i − 1) th solution using Eq. ( 23).
Update the CD of the (i +1) th solution using Eq. ( 24).End For Within the domain of optimizing solution diversity in MOO, various strategies have been explored to enhance the distribution of solutions across the Pareto Front.Among these, the authors in Ref. [89] proposed the Crowding Distance Elimination (CDE) method, specifically designed to improve solution diversity.Concurrently, another set of researchers in Ref. [88] introduced the Dynamic CD-based Diversity Maintenance Approach (DCD-DMA), with a particular focus on expanding solution diversity more uniformly across the objective space.A key feature shared by CDE, DEBCD, and DCD-DMA is their strategic focus on solutions exhibiting the minimum CD.The removal of such solutions triggers a re-evaluation and potentially a reordering of the remaining solutions based on their CD.The primary difference among these methodologies lies in their approach to recalibrating CDs post-removal of a solution.minimizes computational overhead but also enhances the method's efficiency, particularly in scenarios involving extensive solution sets.This efficiency proves to be especially advantageous in maintaining diversity within large solution sets, streamlining the dynamic management of solution distribution across the PF.By adopting DEBCD, practitioners can ensure a more computationally efficient and effective approach to preserving solution diversity in MOOPs, thereby facilitating a more balanced exploration and exploitation of the search space.

Formulation of the proposed algorithm
The methodology of MOO initiates with the creation of an initial population.During this phase, the values of the objective functions are evaluated, and solutions that are non-dominated according to NDS criteria are identified.Following this, the equations governing the RCOA, specifically Eq. (1) through 6, are applied to derive a new set of population members.The selection process for candidates advancing to the subsequent generation involves an analytical comparison of their performance across both previous and current evaluations.In scenarios involving two distinct populations labelled P 1 and P 2 , the determination of dominance is conducted through a comparative analysis of their respective objective function outcomes, denoted as f 1 and f 2 .The dominance of P 1 over P 2 is established if every component of f 1 is either equivalent to or better than (lower in cases where minimization is the goal) its corresponding component in f 2 , with at least one element of f 1 being strictly superior (lower).Solutions that any other candidate does not outperform are classified as NDS solutions and are accordingly assigned a primary rank.Subsequent to the extraction of these top-ranked nondominated solutions, the process iteratively continues to unveil the next tiers of non-dominated solutions, assigning them increasing ranks until every agent within the population is ranked.
The integration and ranking of both existing and newly identified agents are based on their dominance ranks, prioritizing those with lesser degrees of domination.The selection procedure for the ensuing generation commences by accumulating the highest-ranked non-dominated solutions into a new, initially empty collection.Should the aggregate of agents within this collection exceed the designated population limit, the DEBCD method is employed to restrict the agent selectively count within the new collection.The DEBCD strategy is instrumental in managing the set of NDS solutions to adhere to the specified population quota while simultaneously promoting diversity, characterized by a broad and uniform distribution of solutions.Consequently, n solutions exhibiting the most significant CD values are retained for inclusion in the next generation.In instances where the count of rank-1 NDS solutions falls short of the population threshold, solutions with rank-2 designation are incorporated into the collection, a process that is iteratively repeated until the newly formed set aligns with the population size criteria.This NDS-centric selection mechanism is visually represented in Fig. 4, for clarity.Further, Algorithm 3 delineates the pseudocode for implementing the MORCOA strategy.The procedure for MORCOA is initiated with predetermined parameters such as termination criteria Maxit and population size N.An initial random population P 0 within feasible limits is established and subjected to assessment.Subsequently, DEBCD and NDS methodologies are applied to P 0 , leading to the formation of a new population P j , which is then amalgamated with P 0 to generate population P i .This combined population P i is organized based on DEBCD and NDS rankings, from which the optimal solutions are selected to promote a new parent population.This cyclical process is reiterated until the specified termination criteria are satisfied, ensuring a systematic and efficient progression towards the optimization objectives.The operational flow of MORCOA is illustrated in Fig. 5 to provide a ( .sin(10πf1) . ( ) 2 ) S. Ravichandran et al. combined population P i and its reassessment also requires O(N ×dim) space.The storage of the population and objective function values dominates the overall space complexity of MORCOA.Therefore, the overall space complexity is Simplifying, the overall space complexity can be expressed as O(N × (dim + M)).

Best Compromise Solutions
The Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) is a multi-criteria decision-making method that identifies solutions that have the shortest distance from the positive ideal solution (PIS) and the farthest distance from the negative ideal solution (NIS) [90,91].When applied to the MORCOA, TOPSIS can be utilized to determine the Best Compromise Solutions (BCS) among the set of Pareto optimal solutions generated by MORCOA.This approach effectively balances between multiple conflicting objectives by selecting solutions that exhibit an optimal compromise.The steps to implement TOPSIS for MORCOA are as follows.
Step 1. Normalization of the Decision Matrix -The first step involves normalizing the decision matrix to eliminate the effects of differing units and scales among the objectives.The normalized value r ij is calculated using Eq.(25).
where x ij is the value of the j th criterion for the i th alternative, and n is the number of alternatives.
Step 2. : Determination of the PIS and NIS -The PIS A + is the solution with the best value for each criterion among the set of Pareto optimal solutions, and the NIS A − is the solution with the worst value for each criterion.They are defined as follows.
Step 3. Calculation of the Distance to PIS and NIS -The next step is to calculate the distance of each Pareto optimal solution from the PIS and NIS.The distance from the PIS D + i and the distance from the NIS D − i for each solution i are computed as follows.
where m is the number of criteria, r + j and r − j are the values of criterion j for the PIS and NIS, respectively.
Step 4. Calculation of the Relative Closeness to the Ideal Solution -The relative closeness of each solution to the ideal solution is calculated to identify the best compromise solution.The relative closeness C i of solution i is given by: Solutions with higher C i values are considered better, as they are closer to the PIS and further from the NIS.
Step 5. Selection of the BCS -The solution(s) with the highest relative closeness to the ideal solution C i is selected as the best compromise solution.
By implementing TOPSIS in MORCOA, decision-makers can systematically evaluate and select the most suitable compromise solution from the set of Pareto optimal solutions, effectively addressing the trade-offs between multiple objectives.This method provides a clear, mathematical framework for navigating the complex decision space characteristic of multi-objective optimization problems.

Results and discussions
This section explains the research findings and the experimental analyses undertaken to evaluate the MORCOA's efficacy and demonstrate its capabilities.To ensure the robustness and statistical validation of the outcomes, each experimental run was independently conducted 30 times, with the experiments configured to a population size of 40 and a cap of 1000 iterations [92].The computational work was carried out on MATLAB R2023b, hosted on a system powered by an Intel Core i5 processor, clocked at 2.5 GHz, and supported by 16 GB of RAM.The assessment of MORCOA's performance was systematically approached through the employment of three distinct categories of benchmark test functions.
(i) Zitzler-Deb-Thiele (ZDT) Benchmark Suite [41]: This includes a series of test functions, specifically ZDT1 through ZDT4 and ZDT6.These functions are detailed in Table 1 and are characterized by having dual objective functions, mirroring the dual-objective nature prevalent in Pareto optimization, especially within engineering domains.(ii) Deb-Thiele-Laumanns-Zitzler (DTLZ) Benchmark Functions [93]: Distinguished by their scalability, the DTLZ suite permits the extension to an arbitrary number of objectives, rendering it an indispensable asset for exploring challenges in "many-objective" optimization.These functions are detailed in Table 2. (iii) Real-World Engineering Design Problems [94,95]: This category comprises five practical engineering challenges, namely the 4-bar truss design (RW1), disk brake design (RW2), speed reducer design (RW3), welded beam design (RW4), and cantilever beam design (RW5).These problems provide a concrete framework for appraising MORCOA's practical utility and performance in real-world engineering contexts.
The design of honeycomb heat sinks (RW6) is particularly highlighted as a distinct subsection within this study [96].This specific problem encapsulates the intricacies and challenges associated with optimizing thermal management solutions in engineering.Through the honeycomb heat sink design problem, the study showcases MORCOA's adeptness in navigating complex design parameters to achieve optimal thermal performance, underlining the algorithm's relevance and adaptability to engineering optimization tasks.

Results -ZDT and DTLZ functions
This research undertook a comprehensive evaluation of the effectiveness of various optimization algorithms, utilizing a suite of performance indicators to ensure a thorough analysis.These indicators include HV, GD, IGD, SD, SP, and RT.For a more intuitive hold of these metrics, Fig. 2 provides detailed explanations and visual representations of each.In the proposed study, all algorithms are executed individually 30 times, enabling a fair statistical evaluation.Wilcoxon signed-rank test (WSRT) was conducted at a significance level of 0.05 on which the respective algorithms have a better, a worse and an equal performance of MOIPSO, NSGA-II, MOGBO, and MOAGDE with respect to the proposed MORCOA.The outcomes were denoted using the symbols "+/− / = ," representing superior, inferior, and equivalent performance, respectively.This statistical method was applied to contrast the performance of various optimizers across a range of test scenarios, aiming to identify conditions under which these optimizers demonstrated superior, inferior, or comparable effectiveness in solving MOO problems.The comparative study encompassed the proposed MORCOA and several other notable multi-objective optimization algorithms, including (i) Multi-Objective Adaptive Guided Differential Evaluation (MOAGDE): An advanced algorithm that adapts and guides the search process to navigate the solution space efficiently.The MOAGDE locates Pareto optimal solutions for MOOPs with different decision/objective spaces.The MOAGDE was advanced by reshaping the adaptive guided DE algorithm for MOO; (ii) Improved Multi-Objective Particle Swarm Optimization (MOIPSO): A variant of the PSO framework, enhanced for better performance in multi-objective contexts.In each cycle of its operation, the algorithm conducts an assessment of the external population's variety.Based on the outcome, it decides on executing mutation processes on the external population.Furthermore, it selects varying strategies for updating the particle population, guided by the assessed values.This approach is seamlessly integrated with the principles of CD and non-dominant sorting, ensuring a synergistic interaction that enhances the effectiveness of the algorithm; (iii) Multi-Objective Gradient-Based Optimizer (MOGBO): Utilizes Newton's search rule to improve the search process, aiming for a more directed approach to finding optimal solutions; (iv) NSGA-II (Non-dominated Sorting Genetic Algorithm II): NSGA-II stands out in the domain of MOO for its rapid and superior strategy.It utilizes a non-dominated sorting method to rank solutions according to their dominance hierarchy.Esteemed for its capability to rapidly approach the Pareto front, NSGA-II simultaneously ensures the preservation of solution diversity, marking its efficiency and effectiveness in optimization tasks.This research meticulously compares the performance of MORCOA against these algorithms, providing a comprehensive overview and  The summarized performance of various algorithms on ZDT and DTLZ test functions reveals MORCOA as a standout performer.In the context of GD metrics, as listed in Table 3, MORCOA demonstrates superior results across the board, showcasing the smallest values in the majority of the test cases, indicative of its effectiveness in minimizing the generational distance.This is especially evident in the ZDT1 and ZDT2 problems, where MORCOA achieves significantly lower GD values compared to other algorithms like NSGA-II, MOIPSO, MOAGDE, and MOGBO.Its performance suggests an exceptional capability in converging towards the Pareto front while ensuring a diverse solution set.
The insights from the IGD metrics, as listed in Table 4, across ZDT and DTLZ test functions highlight MORCOA's commendable performance.The proposed MORCOA demonstrates a mix of superior, equivalent, and, in fewer cases, inferior performances compared to other algorithms like NSGA-II, MOIPSO, MOAGDE, and MOGBO.Specifically, MORCOA shows outstanding results in ZDT2 and ZDT6 problems, where it outperforms its counterparts with the lowest IGD values, signalling its strength in maintaining a good balance between convergence and diversity towards the Pareto front.
The HV metric analysis, as listed in Table 5, across ZDT and DTLZ test functions indicates that MORCOA exhibits commendable performance.While it does not universally outperform all other algorithms in every test, its results are consistently competitive, often    The analysis of SD metrics, as listed in Table 6, for ZDT and DTLZ test functions highlights MORCOA's commendable performance in maintaining a balanced distribution of solutions across the objective space.While the performance varies across different problems, MORCOA consistently shows either superior or competitive results, particularly shining in scenarios where maintaining a diverse set of solutions is crucial.For instance, in ZDT2 and ZDT4 problems, MORCOA demonstrates significantly better spread metrics compared to other algorithms like NSGA-II, MOIPSO, MOAGDE, and MOGBO, indicating its effectiveness in dispersing solutions evenly along the Pareto front.This is further evidenced in the DTLZ problems, where MORCOA again performs well, especially in DTLZ2 and DTLZ3, showcasing its adaptability and robustness in handling complex multi-objective optimization tasks.
The analysis of SP metrics, as listed in Table 7, across ZDT and DTLZ test functions, reveals an improved performance profile for MORCOA, highlighting its strengths and areas of competitiveness.In some instances, such as ZDT2 and ZDT6, MORCOA showcases superior performance, achieving the lowest spacing metrics, which indicate a more uniform distribution of solutions across the objective space.This is particularly significant in optimization problems where a well-distributed set of solutions is essential for effective decision-making.However, the performance of MORCOA is not uniformly dominant across all test functions.In scenarios like ZDT1, ZDT3, and DTLZ7, its spacing metrics suggest areas where its performance is on par with or slightly lagging behind other algorithms.This mixed outcome underscores the complexity of multi-objective optimization tasks and the challenge of achieving optimal performance across diverse problem sets.MORCOA's performance in spacing metrics, while varied, demonstrates its ability to maintain a competitive edge in ensuring uniform solution distributions in several key test scenarios.
The RT metrics, as listed in Table 8, from ZDT and DTLZ test functions indicate that MORCOA showcases a strong performance in terms of computational efficiency.Across the various problems, MORCOA consistently demonstrates lower run times compared to other algorithms like NSGA-II, MOIPSO, MOAGDE, and MOGBO.Especially notable is its performance in the ZDT1, ZDT2, ZDT3, ZDT4, and ZDT6 problems, where MORCOA not only outperforms its competitors in terms of speed but does so with a considerable margin.In the more complex DTLZ test functions, MORCOA's run times remain competitive, often resulting in faster execution compared to NSGA-II and significantly outperforming MOAGDE.The algorithm shows its ability to maintain efficiency even as the problem dimensionality increases, a crucial factor for practical applications of multi-objective optimization.This consistent pattern of superior run-time performance underscores MORCOA's effectiveness as a computationally efficient algorithm capable of delivering quick solutions without sacrificing quality.
MORCOA emerges as a standout algorithm in the realm of multi-objective optimization based on a thorough evaluation of its performance across a variety of metrics.Its ability is consistently demonstrated through superior results, showcasing its robustness and effectiveness in generating high-quality solutions.Particularly notable is MORCOA's performance in GD and IGD metrics, where it excels in maintaining a balance between convergence towards the Pareto Front (PF) and diversity among solutions.This balance is crucial for effective MOO and is further evidenced in its outstanding results in HV metrics, indicating successful navigation of competing objectives.MORCOA's capabilities extend to SP and SD metrics, where it often outperforms competitors by ensuring a comprehensive exploration of the solution space with evenly distributed solutions.Moreover, its efficiency in the RT metric is increased, demonstrating a swift processing ability that makes it particularly suited for time-sensitive applications.This efficiency, coupled with its adaptability and versatility across various problem types, underscores MORCOA's value in real-world scenarios where   In this research, the newly introduced MORCOA was examined for its efficacy in constructing Pareto fronts across various ZDT benchmark problems.The comparative analysis positioned MORCOA in a leading stance against a suite of established algorithms, including the NSGA-II, MOIPSO, MOGBO, and MOAGDE.Fig. 6, integral to this study, visually delineates the PFs as generated by each algorithm under consideration.Notably, Fig. 6 comprises individual PF plots for each algorithm alongside a collective PF plot, facilitating a direct comparison of their performance metrics.A significant highlight from the analysis was MORCOA's noticeable accuracy in navigating towards the true PF, a superiority that was markedly observable within the more intricate regions of the objective space.Interestingly, this precise convergence was achieved without the necessity for any parameter tuning within the proposed algorithm, underscoring its inherent robustness and adaptability.Moreover, the PF realized through MORCOA was characterized by a more uniform distribution of solutions.This uniformity is indicative of the algorithm's comprehensive exploration capabilities across the objective space, ensuring that a broad spectrum of viable solutions is considered.In contrast, the comparative algorithms frequently exhibited tendencies towards solution clustering or the presence of significant gaps within their generated PFs.Such discrepancies highlight the limitations of these algorithms in adequately spanning the objective space, thereby highlighting the enhanced performance and efficiency of MORCOA in multi-objective optimization scenarios.
In this study, the MORCOA showcased exceptional adaptability and performance across a range of benchmark problems, notably the ZDT and DTLZ test suites, which are renowned for their complexity and varied PF geometries, including convex, concave, and discontinuous shapes.Unlike its counterparts -MOIPSO, MOGBO, NSGA-II, and MOAGDE -MORCOA demonstrated unparalleled consistency in performance across multiple runs, showcasing less variance in solution quality and indicating higher reliability for multi-objective optimization tasks.
MORCOA's adaptability was particularly highlighted in its handling of the DTLZ problems set with tri-objective challenges, where it not only achieved remarkable precision in converging towards the true PF but also excelled in managing the intricate interdependencies between objective functions.This level of precision in convergence, particularly under the complex conditions presented by tri-objective scenarios, sets MORCOA apart from competing algorithms.The algorithm's sophisticated structure was adept at navigating these complexities, leading to Pareto fronts characterized by superior solution distribution and diversity.Unlike the other algorithms, MORCOA produced PFs that were more uniformly distributed, ensuring a comprehensive exploration of the objective space.This uniformity is vital for capturing the entire spectrum of possible trade-offs between competing objectives.Fig. 7 in the study highlights MORCOA's superiority, showing its ability to generate more accurate and diverse PFs compared to its counterparts.This capability is crucial for understanding the full range of solutions available in MOO scenarios.Moreover, despite its progressive optimization features, MORCOA maintained computational efficiency on par with or exceeding that of the compared algorithms.This efficiency is particularly important in scenarios with limited computational resources, highlighting MORCOA's practicality in resource-constrained environments.Overall, MORCOA's robustness, adaptability, and computational efficiency underscore its effectiveness and potential as a powerful tool for tackling complex multi-objective optimization problems.
The visualization of box plots, as shown in Figs. 8 and 9, for various metrics such as RT, SP, SD, HV, GD, and IGD across the ZDT and DTLZ problem sets, respectively, serves as a comprehensive statistical tool to evaluate and compare the performance of the proposed MORCOA against other algorithms.These plots provide a multifaceted view of MORCOA's reliability and efficacy in MOO by     illustrating the range, median, and variability of each metric across all tested problems.The RT box plots depict computational efficiency of MORCOA relative to other algorithms.A lower median and tighter interquartile range (IQR) suggest MORCOA is not only fast but also consistent in its execution time across multiple runs, indicating efficient algorithmic complexity and implementation.SP metrics assess the uniformity in the distribution of solutions along the Pareto front.Box plots showcasing lower medians and smaller IQRs for MORCOA indicate its ability to generate evenly spaced solutions, contributing to its reliability in providing a diverse set of solutions that are well-distributed across the objective space.The SD box plots evaluate the extent to which the solutions cover the true Pareto front.MORCOA's performance in this metric, indicated by lower variability and a median close to optimal values, demonstrates its effectiveness in exploring the objective space thoroughly and capturing a wide range of optimal solutions.HV represents the size of the space covered by the Pareto front relative to a reference point.
Box plots with higher medians and narrow IQR for MORCOA signify its superior ability to achieve solutions that not only converge to the true Pareto front but also cover a substantial portion of the objective space, highlighting its robustness and reliability in capturing the true trade-offs between objectives.GD and IGD box plots provide insights into the convergence and diversity of the solutions, respectively.For MORCOA, lower GD and IGD values indicate a strong convergence towards the true Pareto front while maintaining a diverse set of solutions.This balance between convergence and diversity underscores MORCOA's reliability in accurately and comprehensively addressing MOO problems.The box plot visualization collectively argues for MORCOA's reliability in multi-objective optimization.The algorithm demonstrates not only efficiency and speed (as seen in RT plots) but also excels in providing high-quality solutions that are both diverse and well-distributed (SP, SD), cover a significant portion of the objective space (HV), and closely approximate the true Pareto front with a balanced set of solutions (GD, IGD).The consistency in MORCOA's performance across a wide range of problems, as indicated by the limited variability in the box plots, further cements its reliability and effectiveness as a robust tool for solving complex MOO problems.

Real-world engineering problems
In this comprehensive evaluation, the capabilities of the proposed MORCOA are rigorously examined across a range of difficult  engineering design challenges.These real-world scenarios, representative of the multimodal and multi-objective problems frequently encountered within the engineering domain, serve as a testament to the complexity and practical significance of the tasks at hand.The specific engineering problems scrutinized in this study include (i) 4-bar truss design problem, (ii) welded beam design problem, (iii) disk brake design problem, (iv) cantilever beam design problem, (v) speed reducer design problem, and (vi) honeycomb heatsink design problem.Each of these design problems showcases a distinct aspect of engineering optimization, pushing the boundaries of MORCOA to balance multiple objectives within a constrained solution space effectively.The algorithm's performance across these diverse scenarios reveals its utility and effectiveness as a powerful tool for attempting the multi-dimensional optimization tasks inherent in engineering applications.Through this rigorous assessment, MORCOA's proficiency in navigating the complex landscape of engineering design problems is intensely demonstrated, reinforcing its value as an adaptable and reliable benefit for multi-objective optimization problems.

Four-bar truss design problem
In the analysis of the 4-bar truss optimization scenario delineated in source [73], the objective functions f 1 and f 2 target the minimization of structural volume and nodal displacement, respectively.This optimization conundrum is critical within the domain of structural engineering, emphasizing the requisites for material utilization efficiency and structural stability.The problem is defined by four design variables x 1 to x 4 , corresponding to the cross-sectional areas of truss members 1 through 4, which are integral in modulating both the aggregate structural volume and displacement under operational loads.The minimization of the structural volume, denoted by f 1 , entails reducing the material footprint by optimizing the variables related to the cross-sectional areas and the linear dimensions of the truss members, thus directly affecting the volume, as articulated in Eq. (28).Concurrently, the reduction of nodal displacement, encapsulated in f 2 , is pivotal for ensuring structural efficacy and resilience, also formalized in Eq. ( 28).This dual-faceted optimization strategy underscores the intricate interplay between achieving material efficiency and ensuring structural integrity within engineering design parameters.
The displacement experienced by the truss arises from its structural reaction to applied loads, with the cross-sectional areas of its members playing a significant role in determining this response.This displacement can be quantified through the application of structural analysis equations that take into account the truss's geometry, the properties of the materials used, and the specific loading conditions it faces.The design variables x 1 , x 2 , x 3 , and x 4 represent the cross-sectional areas of the truss's four members, respectively.The values assigned to these variables are critical, as they not only dictate the overall structural volume of the truss but also have a direct impact on its displacement under load.The process of designing the truss is governed by a set of constraints, as outlined in Eq. (29).
These constraints typically encompass limitations on the stress experienced by each member, a cap on the allowable displacement of the truss as a whole, and additional considerations for the prevention of buckling in members under compression.Such constraints ensure that the final truss design maintains both structural integrity and functionality under the prescribed loading conditions.Fig. 10 presents a diagrammatic view of the truss, detailing the layout of its components and delineating how the design variables x 1 , x 2 , x 3 , and x 4 relate to each member.

Disk brake design problem
The disk brake design challenge, characterized by mixed constraints, was delineated by Osyczka and Kundu [97] within the domain of automotive engineering, aiming at the dual optimization of disk brake performance and efficiency, as depicted in Fig. 11.This issue centralizes on reducing the stopping time f 1 and the mass f 2 of the disk brake, essential for augmenting both the safety and operational efficacy of braking systems.The disk's inner radius x 1 plays a pivotal role in determining thermal capacity and structural resilience, while its outer radius x 2 is key to modulating braking torque and heat dissipation rates.Additionally, the engaging force x 3 is intrinsically linked to the exerted braking force, thereby influencing the stopping time.The number of friction surfaces x 4 directly impacts the brake's frictional interface, affecting both stopping time and thermal efficiency.
Reducing the stopping time is crucial for vehicular safety; facilitating an efficient halt of the vehicle over minimal distances while decreasing the brake mass contributes to the overall reduction in vehicle weight.This reduction is beneficial for enhancing fuel economy and vehicle performance.The defined objective functions and associated constraints are outlined as follows.
Minimize : Subjected to: The limits are provided in Eq. ( 32). 55 The Pareto optimal front for this problem is convex and continuous, showing a direct trade-off between the stopping time and the brake mass.This means that as one objective gets better, the other worsens consistently and measurably.This issue serves as a typical example of MOO in engineering, highlighting the need to find a balance between performance (stopping time) and practical aspects (mass), considering the limits set by physical and mechanical realities.

Speed reducer design problem
The speed reducer design problem, detailed in Ref. [65] and depicted in Fig. 12, represents a quintessential multi-objective optimization issue within mechanical engineering, aiming to reduce both the weight f 1 and stress f 2 of a speed reducer.These Fig. 14.Cantilever beam design.
S. Ravichandran et al. objectives are vital for improving the gear system's efficiency and longevity.The design involves seven decision variables that significantly influence the system's performance and structural resilience.Variables such as the gear face width x 1 and teeth module x 2 dictate the gear's physical dimensions and operational behaviour, including load capacity and gear teeth spacing.The number of teeth on the pinion x 3 , an integer variable, defines the gear ratio and the pinion's mechanical leverage.The distances between bearings x 4 and x 5 determine load distribution and the assembly's mechanical stability, while the diameters of shafts 1 and 2 (x 6 and x 7 ) are critical for ensuring structural integrity and torque transmission capabilities of the respective shafts.This problem exemplifies the balance needed between performance enhancement and practical design constraints in engineering optimizations.
The objective of this problem is to reduce the speed reducer's weight to increase efficiency and lower material costs while also minimizing stress to ensure the gear system's durability and reliability under load.The objective functions and constraints are outlined as follows.
Minimize : Subjected to: ) − 1 The related parameters of this design are given as follows.σ = 504000

Welded beam design problem
The welded beam design problem, introduced by Deb [82], stands as a significant case in structural optimization, focusing on minimizing fabrication cost f 1 and beam deflection f 2 .These objectives are vital for ensuring cost-effectiveness and maintaining structural integrity.Illustrated in Fig. 13, the design involves four decision variables: the weld thickness x 1 , affecting joint strength and cost; the clamped bar length x 2 , influencing beam dimensions and stability; the bar height x 3 , crucial for load-bearing and deflection, and the bar thickness x 4 , impacting structural behaviour.The challenge aims to reduce material and labour costs while ensuring the beam's structural efficiency and durability under operational loads with specified objective functions and constraints, as follows.S. Ravichandran et al.

Cantilever beam design
This study addresses the cantilever beam design problem, a core issue in structural engineering, characterized by three design variables and two objectives under two specific loading scenarios, depicted in Fig. 14 [65].The variables include the beam's width x 1 , height x 2 , and total length x 3 , all measured in inches.The goals are to minimize the beam's weight for material efficiency and application versatility and to reduce the strain energy, which influences the beam's deflection and structural integrity.The loading conditions considered are: (i) an end load causing maximum bending stress σ, and (ii) a torque load inducing shear stress τ at the beam's root.The optimization challenge integrates these dimensions and objectives within the framework of the specified loading conditions, aiming to enhance the beam's performance and durability.The fitness functions and constraints are as follows.
Minimize : Subjected to: The related parameters are given as follows.The maximum allowed stress S y is selected as 300 × 10 3 in kPa, the elastic modulus E is    selected as 207 × 10 6 in kPa, ρ is selected as 7800 in kg/m 3 , the load P is selected as 1 in kN, and the maximum allowed deflection Δ max is selected 0.005 in mm.

Honeycomb heat sink design problem
Honeycomb assemblies have garnered widespread recognition across multiple engineering disciplines for their exceptional attributes and versatile functionality.These assemblies stand out in applications demanding efficient thermal regulation, absorption of impact energy, and support against structural loads.The unique combination of thermal properties and geometric configuration of honeycomb structures yields several key advantages [96].Firstly, their significant surface-to-area ratio amplifies their utility in heat transfer applications, enabling superior efficiency in scenarios where the surface area is paramount.Secondly, the architectural layout of honeycomb structures promotes a streamlined fluid flow, presenting an optimal solution for scenarios requiring the preservation of fluid dynamics with reduced energy dissipation.Lastly, the construction materials, commonly including metals such as aluminium, possess high thermal conductivity, ensuring effective heat dispersion.
Central to heat sink optimization is the principle of thermal resistance, which measures the challenge heat encounters during transfer through materials.Optimizing a heat sink involves minimizing this resistance, thereby enabling quicker heat dissipation.Materials like aluminium are often chosen for their low thermal resistance, balancing efficiency with cost and weight considerations.The unique geometry of honeycomb structures significantly increases the surface area available for convective heat transfer, as indicated by the Nusselt number (Nu).This dimensionless metric assesses the convective heat transfer rate relative to conductive heat transfer, guiding the design process to maximize Nu through strategic modifications of the heat sink's geometry and surface characteristics.Fluid dynamics principles, particularly the friction factor (ff), are crucial in evaluating the resistance fluid experiences flowing through the heat sink.Optimizing design to reduce the friction factor ensures efficient fluid flow with minimal energy loss, which is crucial for systems where air or liquid serves as the cooling medium.The choice of material impacts not only the thermal conductivity of the heat sink but also its durability, weight, and overall cost.Aluminium is frequently selected for its excellent thermal conductivity and cost-effectiveness, highlighting the need for materials that enhance thermal performance without compromising structural integrity or economic viability.Adjusting the honeycomb structure's geometry, including cell size, shape, and arrangement, directly affects thermal performance.Design variables such as cross-sectional areas and dimensions of components are meticulously optimized to balance thermal efficiency with physical and mechanical constraints.Design optimization also accounts for the specific heat load and ambient conditions expected in operation.This ensures that the heat sink design is robust and capable of performing under varied and potentially extreme conditions.
In the background of heat sink and heat exchanger design, honeycomb structures, as elucidated in sources [92], demonstrate exceptional proficiency.Their structural design facilitates stellar heat transfer performance while simultaneously optimizing space utilization and curbing energy wastage.When crafting an exemplary heat sink, especially one adorned with hexagonal aluminium fins, the analysis of five decision parameters becomes imperative.These parameters, encompassing fin dimensions, inter-fin spacing, material gauge, and fin elevation, critically impact the heat sink's thermal and frictional performance.Evaluation of such a heat sink's efficiency hinges on analyzing two essential objective functions: the friction factor (ff) and the Nusselt number (Nu).The Nu, a dimensionless metric, gauges the effectiveness of convective over conductive heat transfer at the fluid boundary, with higher values  indicating enhanced convective heat transfer efficiency.Conversely, the ff metric quantifies the frictional resistance encountered by the fluid navigating through the heat sink, underscoring the importance of minimizing friction to bolster fluid flow efficiency and minimize energy expenditure.Through meticulous optimization of these parameters, the aim is to engineer a heat sink design that optimally balances heat transfer enhancement against frictional resistance, thereby elevating the efficacy and efficiency of thermal management systems.The objective functions are detailed as follows.

Results of real-world engineering problems
Table 9 within this study presents an in-depth statistical analysis comparing the MORCOA against a range of alternative optimization algorithms, with a particular focus on their efficacy in solving diverse engineering design problems.The primary metric employed for this evaluation is the GD, a measure that assesses the proximity of an algorithm's solutions to the true Pareto front, thereby indicating its convergence capability.The data compiled in Table 9 clearly demonstrates MORCOA's superior performance relative to the other considered algorithms across a variety of engineering design scenarios.This superiority is quantitatively underscored by MORCOA's consistently lower GD values across the board.Such lower GD values are indicative of MORCOA's superior convergence properties, showcasing its ability to closely approximate the true Pareto front more effectively than its counterparts in these complex engineering optimization tasks.In addition, Friedman's ranking test (FRT) values are also provided.
In this section, Table 9 provides a detailed statistical benchmark, offering a comparative analysis of MORCOA and other optimization algorithms in tackling various engineering design challenges.The primary metric for this evaluation is the GD, which measures how closely an algorithm's generated solutions approach the true Pareto front, indicating the algorithm's convergence efficiency.The data in Table 9 shows that MORCOA achieves consistently lower GD values, highlighting its superior convergence capabilities and proficiency in more accurately approximating the true Pareto front within complex engineering optimization landscapes.Additionally, Fig. 18.RT analysis of all algorithms for real-world problems.

S. Ravichandran et al.
Table 9 evaluates the algorithms using the IGD metric, which assesses how effectively these algorithms solve a series of engineering design problems.The analysis reveals that MORCOA consistently outperforms other prominent algorithms such as NSGA-II, MOIPSO, MOGBO, and MOAGDE across five distinct real-world engineering challenges.This performance level underscores MORCOA's competitive advantage and robustness, especially in contexts where performance margins among different algorithms are relatively narrow.Table 9 further broadens the scope of our analysis by incorporating additional pivotal performance indicators: SD, SP, and HV, each offering unique insights into the efficacy of the optimization algorithms under consideration: The comparative results reveal a notable trend: MORCOA, the algorithm introduced in this study, exhibits superior performance in the majority of cases when evaluated against the backdrop of these metrics.Specifically, MORCOA outperforms its competitors in five out of six scenarios for the Spread metric.This consistent superiority underlines MORCOA's exceptional capability not only to discover a broad spectrum of optimal solutions but also to ensure that these solutions are evenly distributed across the entirety of the Pareto front.MORCOA also shows outstanding performance in the Spacing metric, ensuring an equitable distribution of solutions across the Pareto front, which is critical in multi-objective optimization to avoid solution clustering and guarantee comprehensive exploration of the solution space.In most of the scenarios assessed, MORCOA surpasses its competitors, further highlighting its efficiency in navigating multi-objective optimization challenges.The extensive statistical analysis provided across these tables illuminates the durability and adaptability of MORCOA in managing a variety of complex engineering design problems, demonstrating superior performance over other established algorithms across different evaluative metrics.This detailed analysis showcases MORCOA's dominance in striking an optimal balance between finding a comprehensive set of solutions and ensuring their equitable distribution, thereby enhancing the overall decision-making process in multi-objective optimization scenarios.
To definitively demonstrate that the Pareto fronts generated by MORCOA are more aligned with the optimal solutions than those produced by alternative algorithms, this study includes detailed visual evidence in Fig. 16.This illustration compares the true Pareto fronts with those derived through various algorithms for the engineering design challenges addressed.A meticulous examination of these visuals reveals several critical insights into MORCOA's performance efficacy.Specifically, Fig. 16 clearly illustrates that the Pareto fronts realized by MORCOA more closely approximate the true Pareto fronts than those generated by other algorithms, signalling a higher precision in MORCOA's solution set.This proximity to the true Pareto fronts evidences MORCOA's superior capability in identifying solutions that adeptly balance the inherent trade-offs among the multiple objectives characteristic of these engineering tasks.
MORCOA's ability to produce well-distributed solution sets across the objective space stands out distinctly.Contrary to other algorithms, which may exhibit solution clustering or significant gaps, MORCOA spans a more extensive portion of the objective space, reflecting a deep and varied exploration of possible solutions.This broad distribution is imperative for encapsulating a comprehensive array of trade-offs, thus offering decision-makers a wide selection of optimal choices.The consistent achievement of MORCOA across diverse engineering design challenges, as depicted in Fig. 16, affirms its robustness and versatility.Regardless of the complexity of the structural or mechanical system in question, MORCOA invariably yields Pareto fronts that are markedly congruent with the true fronts.The graphical presentations in Fig. 16 provide a vivid and intuitive validation of MORCOA's operational effectiveness.
In this comprehensive study, the inclusion of box plots, as shown in Fig. 17, for all metrics across the entire set of problems provides a detailed statistical overview, further substantiating the reliability and competency of the MORCOA when compared with state-of-theart algorithms.These box plots serve as a visual testament to MORCOA's performance, offering insights into the central tendency, dispersion, and outliers of the results obtained for each metric, thereby facilitating a granular analysis of its effectiveness across various engineering design challenges.The box plots meticulously detail the performance distribution for metrics such as GD, IGD, Spread, Spacing, and HV.Through these plots, MORCOA's ability to consistently achieve lower GD and IGD values is visually highlighted, demonstrating its superior convergence towards the true Pareto fronts.Similarly, the Spread and Spacing metrics depicted in the box plots reveal MORCOA's effectiveness in ensuring a diverse yet evenly distributed set of solutions across the Pareto front, an essential aspect of multi-objective optimization that underscores the algorithm's ability to explore the solution space comprehensively.Moreover, the HV metric, as visualized through the box plots, underscores MORCOA's capacity to cover a significant volume of the objective space, indicating not only the quality but also the quantity of the solutions it generates.This extensive coverage affirms MORCOA's adeptness in identifying a broad spectrum of optimal solutions, thereby providing decision-makers with a wide range of viable options to consider.
Runtime (RT) analysis is crucial for comparing the efficiency of the proposed algorithm, assessing its scalability with problem complexity, and ensuring practical applicability in real-time scenarios.To better understand the performance of all algorithms, Table 10 and Fig. 18 are provided to visualize their runtimes across the different real-world problems.
It is observed from Fig. 18 and Table 10 that MORCOA generally exhibits the lowest runtime across most problems, and it shows significant efficiency, particularly in RW1 and RW2.NSGA-II consistently has higher runtimes compared to other algorithms and demonstrates less efficiency, particularly in RW2 and RW3.MOAGDE's runtimes are generally in the middle range, higher than MORCOA but lower than NSGA-II for most problems and exhibit balanced performance across all problems.MOIPSO exhibits high runtimes, especially in RW3, RW4, RW5, and RW6 and indicates less efficiency in handling these specific problems.MOGBO is similar to MOIPSO, shows higher runtimes in later problems and demonstrates a trend of increased runtime in more complex problems.The RT analysis highlights the efficiency of MORCOA in comparison to other algorithms.MORCOA consistently shows lower runtimes, indicating its suitability for real-world problems where computational efficiency is crucial.
The aggregate analysis provided by box plots and RT analysis, across all considered metrics and problems, provides a depiction of MORCOA's reliability and superior performance.Compared to other contemporary algorithms, MORCOA demonstrates enhanced efficiency, robustness, and adaptability in addressing the multifaceted requirements of complex engineering design problems.The comprehensive visual and statistical evidence laid out in the box plots and RT analysis clearly supports the statement that MORCOA S. Ravichandran et al. stands as a more competent solution in the landscape of multi-objective optimization algorithms, offering a reliable and effective tool for attempting the complicated trade-offs and challenges inherent in engineering optimization tasks.

Limitations of this study
In Tables 6 and 7, the performance of MORCOA in terms of Spacing and Spread metrics shows some variability, especially when compared to state-of-the-art algorithms.These limitations can be attributed to several factors.
• While effective in maintaining diversity, DEBCD may sometimes result in suboptimal spacing.This occurs because the elimination process might remove solutions that could have contributed to a more even distribution.• The ZDT and DTLZ test suites encompass a wide range of problem complexities, including deceptive Pareto fronts and varying degrees of difficulty.MORCOA's transient response behaviour is highly effective for certain problem types but may struggle with others, leading to variations in spacing and spread performance.• Achieving this balance perfectly across all test problems is challenging.In some cases, the algorithm might over-exploit certain regions, leading to clustered solutions and poor spacing.• Although MORCOA is designed to handle a range of optimization problems, its efficiency and effectiveness may diminish with increasing problem scale and dimensionality.Large-scale problems require more computational resources and sophisticated handling of multi-dimensional spaces, where MORCOA might struggle to maintain its performance.
To enhance the robustness and efficiency of MORCOA, the following strategies could be implemented: (i) Introducing mechanisms that dynamically adjust the resistance and capacitance parameters during the optimization process based on real-time feedback can improve the algorithm's adaptability to different problem complexities; (ii) Incorporating advanced diversity preservation techniques, such as clustering-based methods or niche-preserving strategies, can help maintain a well-distributed Pareto front and improve spacing metrics; (iii) Combining MORCOA with other optimization techniques that excel in different areas of MOOPs, such as differential evolution or particle swarm optimization, could leverage the strengths of multiple algorithms, leading to better overall performance; (iv) Developing more efficient algorithms and parallel processing techniques can help MORCOA scale better with increasing problem size and complexity, making it more applicable to large-scale optimization problems.
While MORCOA shows promise in addressing multi-objective optimization challenges, acknowledging and addressing its limitations is crucial for further development.By implementing adaptive parameter tuning, enhancing diversity preservation, exploring hybrid approaches, and improving scalability, MORCOA's performance can be significantly enhanced, making it a more robust and versatile tool for complex optimization tasks.

Conclusions and future scopes
This investigation has rigorously validated the efficacy of the proposed MORCOA through its application to a series of numerical benchmarks and real-world engineering challenges, including the complex honeycomb heat sink design problems.Across these various testing scenarios, MORCOA has consistently demonstrated a superior performance profile when benchmarked against state-of-the-art MOO algorithms.Key findings from this study include MORCOA's exceptional ability to converge towards the true Pareto fronts, significantly outperforming competitors such as NSGA-II, MOIPSO, MOGBO, and MOAGDE.This was particularly evident in engineering design problems where MORCOA not only identified different optimal solutions but also ensured these solutions were evenly distributed across the Pareto front, showcasing its robust capability in achieving a comprehensive exploration of the solution space.In the context of the honeycomb heat sink design optimization, MORCOA's superior performance was highlighted by its proficient handling of the design's complex thermal and structural requirements.By effectively balancing objectives like minimizing thermal resistance and optimizing fluid flow, MORCOA illustrated its adaptability and precision in solving the complicated design spaces typical of real-world engineering problems.The statistical analysis, enriched by the use of performance metrics such as GD, IGD, Spread, Spacing, and HV, provided a quantitative foundation for MORCOA's evaluation.The algorithm's dominance in these metrics across various problem sets reinforces its status as a highly reliable and efficient tool for multi-objective optimization tasks.The study's comprehensive approach, combining theoretical insights with empirical analysis, has clearly established MORCOA as a valuable asset in the field of optimization.Its proven capability to deliver high-quality, diverse, and evenly distributed solutions across a wide range of optimization problems positions MORCOA as a leading choice for engineers and researchers facing complex multi-objective design and optimization challenges.
Given the promising results obtained in this study, several avenues for future research emerge: (i) While MORCOA has demonstrated robust performance, exploring potential algorithmic enhancements could further improve its efficiency.This might include integrating adaptive mechanisms for parameter tuning or exploring hybridization with other optimization techniques to enhance its exploratory and exploitative capabilities; (ii) Extending the application of MORCOA to other complex, real-world problems beyond engineering design could provide insights into its adaptability and efficiency in different domains, such as biomedical, financial, and environmental engineering; (iii) Investigating the scalability of MORCOA, particularly in handling problems with a higher number of objectives and decision variables, would be valuable; (iv) As new optimization algorithms are developed, conducting comparative studies remains crucial.Future research should continue to benchmark MORCOA against emerging algorithms to ensure its competitiveness and identify areas for improvement; (v) Developing decision support systems incorporating MORCOA's optimization capabilities could be beneficial.

Definition 4 .
Pareto Optimal Front[82] -The Pareto Optimal Front, or PF, represents the mapping of the Pareto Optimal Set in the objective space.It is a graphical depiction that illustrates the trade-offs between the different objectives that the Pareto optimal solutions embody.Each point on the Pareto Front corresponds to a Pareto optimal solution, showcasing the best trade-offs achievable among the objectives.The Pareto Front serves as a valuable tool for decision-makers, providing a visual overview of the possible outcomes and facilitating informed decisions based on the relative importance of each objective.The PF is mathematically expressed as follows.

Fig. 3 .
Fig. 3. Modification of CD of solutions nearer to removed one using dynamic-elimination CD for two objectives.(a) CDs of x i− 1 and x i+1 before x i is removed, (b) CDs of x i− 1 and x i+1 after x i is removed.

Algorithm 2 .
DEBCD Method for External Archive Maintenance n: Total number of solutions in the fixed external archive.p: Number of solutions to be removed to maintain the fixed external archive size.For each target solution j from 1 to p: Sort the solutions by their CD for i in the range (2, 3, …, n − 1).

S
.Ravichandran et al.   leading or closely matching the best-performing algorithms.Particularly noteworthy is MORCOA's superior performance in specific problems such as ZDT2, where it achieves the highest HV value, suggesting an effective balance between solution diversity and convergence toward the Pareto front.Similarly, in DTLZ test functions, MORCOA's performance is robust, often matching or closely following the leading scores.

S
.Ravichandran et al.

S
.Ravichandran et al.

S
.Ravichandran et al.

S
.Ravichandran et al.

S.
Ravichandran et al.

•
If the current iteration equals the maximum allowed iterations, conclude the optimization and return the best solution.Otherwise, repeat from step 2.
The pseudocode of the RCOA is provided in Algorithm 1. Algorithm 1. Pseudocode for the RCOA Initialize Parameters: For each individual in the population Xi, where i = 1, 2, …, N and N is the size of the population and set initial values for time t and time constant τ.Initial Fitness Evaluation: Calculate the fitness of each member of the population.Identify transient-state Xts ̅→ and steady-state Xss ̅→ components.While It < Maxit do (current iteration is less than the maximum iterations):

Table 1
understanding of the algorithm's implementation.Execution Procedure for MORCOA Initialize RCOA parameters and the initial population P0(t = 0) within the feasible solution space S. 1. Evaluate the objective function values F for the initially generated population P0. 2. Apply NDS to P0 to determine non-dominated ranks and identify distinct fronts.3.For each identified front, calculate CD values using Eqs.(23) and (24).4. Utilize the RCOA as outlined in Algorithm 1 to update the solutions, resulting in population Pj. 5. Merge the initial population P0 with the updated population Pj to form a combined population Pi = P0 ∪ Pj. 6. Reassess the combined population Pi within the objective space F. 7. Organize Pi according to non-dominated rankings and associated CD values.8. Update the initial population P0 by selecting N top-ranked members from Pi based on their non-dominated rank and CD.Return: Updated population representing the best optimal fronts.Initializing the population P 0 and setting RCOA parameters takes O(N ×dim) time.Evaluating the objective function values for the initial population P 0 requires O(N ×dim) time.Applying NDS to P 0 to determine non-dominated ranks and identify distinct fronts takes O ( M ×N 2 ) time, where \ M is the number of objectives.Calculating the crowding distance values for each identified front using DEBCD takes O ( N 2 ) time.Utilizing RCOA to update the solutions, resulting in population P j , takes O(N ×dim) time per iteration.For Maxit iterations, this takes O(N ×dim ×Maxit) time.Merging the initial population P 0 with the updated population P j to form a combined population P i takes O(N) time.Reassessing the combined population P i within the objective space F requires O(N ×dim) time.Reapplying NDS on P i takes O ( M ×N 2 ) time.Selecting N top-ranked members from P i based on their non-dominated rank, and CD takes O(N log N) time.Considering all steps, the overall time complexity of MORCOA is dominated by the non-dominated sorting and RCOA updates.Hence, the overall time complexity is O(Maxit ×N ×dim) + O List of ZDT benchmark test functions.
S.Ravichandran et al.comprehensive( M ×N 2 ) + O ( N 2 ) and simplified as O ( Maxit × N × dim + M × N 2 ) .3.6.2.Space complexityStoring the population P 0 requires O(N ×dim) space.Storing the objective function values for each individual in the population requires O(N ×M) space.Storing crowding distance values and non-dominated ranks for each individual requires O(N) space.The

Table 2
Lists of three objective DTLZ benchmark test functions.

Table 3
GD metrics of ZDT and DTLZ test functions.

Table 4
IGD metrics of ZDT and DTLZ test functions.

Table 5
HV metrics of ZDT and DTLZ test functions.

Table 6
SD metrics of ZDT and DTLZ test functions.

Table 7
SP metrics of ZDT and DTLZ test functions.

Table 8
RT metrics of ZDT and DTLZ test functions.

Table 9
All performance metrics for RW engineering problems.

Table 10
Average RT values of all algorithms after 30 individual runs.
S.Ravichandran et al.

Table A1
Control parameters of all selected algorithms